Magnon Heat Transport in doped La 2 Cu04 
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We present results of the thermal conductivity of La2Cu04 and Lai.sEuo^CuCU single-crystals 
which represent model systems for the two-dimensional spin- 1/2 Heisenberg antiferromagnet on 
a square lattice. We find large anisotropies of the thermal conductivity which are explained in 
terms of two-dimensional heat conduction by magnons within the Cu02 planes. Non-magnetic Zn 
substituted for Cu gradually suppresses this magnon thermal conductivity K mm . A semiclassical 
analysis of K mag is shown to yield a magnon mean free path which scales linearly with the reciprocal 
concentration of Zn-ions. 
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Apart from their intimate relation to the physics of 
the copper-oxide high-temperature superconductors low 
dimensional quantum magnets are at the center of many 
novel phenomena. Among them are magnetic contribu- 
tions to the thermal conductivity k of an unprecedentedly 
large magnitude, which have been discovered recently in 
the quasi one-dimensional (ID) spin ladder compound 
(Sr,Ca,La)i 4 Cu240 4 i 013- 

Similar findings of large and 
quasi ID magnetic heat conduction have been reported 
also for cuprate spin-chain compounds 0,0- In this con- 
text it has become a key issue whether thermal transport 
due to magnetic degrees of freedom in certain classes of 
spin models may be intrinsically quasi- ballistic, i.e., dissi- 
pationless, thereby leading to high magnetic thermal con- 
ductivities K mag [a, @ ■ Of particular interest is the mean 
free path Z mag of the heat-current carrying magnetic ex- 
citations. Based on a kinetic approach it has been shown 
recently 0, that Z mag for the spin ladder compound is 
of the order of 3000 A . Direct microscopic evaluation of 
the mean free path is an open issue yet, with conflicting 
conclusions published @. 

In this letter we turn to two-dimensional (2D) quan- 
tum magnets by investigating the thermal conductiv- 
ity of La2CuC>4 which realizes the spin- 1/2 Heisenberg- 
antifcrromagnet on the square lattice. The magnetic 
structure consists of Cu 2+ -ions in CuCVplanes extend- 
ing along the crystallographic a- and 6-directions. A 
strong antiferromagnetic intraplanar exchange coupling 
J/Izb ~ 1550 K pi is present whereas the interplanar 
exchange J± is negligible (J±/J~ 10 -5 ) @. 

Our primary focus will be on a broad peak at high 
temperatures (T) in k(T) which we observe if measured 
by a thermal current parallel to the Cu02-planes, but 
which is absent for the perpendicular direction. This is 
consistent with previous results obtained in Ref. j^l ■ Sev- 
eral attempts to explain this peak have been made, e.g. 
scattering processes of acoustic phonons with soft opti- 
cal phonons (JJj or magnons ^lj ■ It has been speculated 
also, that the high-T-maximum could be related to heat 
carried by magnetic excitations However, its origin 



has never been elucidated unambiguously. Here, and em- 
ploying the effects of doping by various types of impuri- 
ties we provide clear evidence in favor of the anisotropic 
heat conduction to be due to 2D magnons indeed. We ex- 
tract the mean free path from our data using an approach 
similar to that for the case of (Sr, Ca, La)i4Cu2404i |]J. 
It will be shown that Z mag does not only scale linearly 
with the reciprocal concentration of magnetic impurities 
which are generated by Zn-ions, but moreover, that it 
is roughly equal to the mean unidirectional distance be- 
tween these ions. 

We have measured k of La2Cu0 4 and Lai.sEuo.2Cu04 
single crystals as well as La2Cui_ z Zn z 04 (z = 0, 0.005, 
0.008, 0.01, 0.02, 0.05) polycrystals as a function of T. 
Stoichiometric oxygen contents were achieved by anneal- 
ing in high vacuum (p < 10~ 4 mbar) for 1-3 hours at 
800° C. The preparation of the samples 01 an d the ex- 
perimental method Q have been described elsewhere. 

In the lower panel of Fig.^we present the thermal con- 
ductivity of single crystalline La2Cu04 if measured par- 
allel (K a b) an d perpendicular (k c ) to the Cu02-planes as 
a function of T. The pronounced anisotropy of k strongly 
resembles that of spin ladder compounds jjj] : k c exhibits 
a T-dependence which is typical for insulating crystalline 
materials with pure phononic thermal conductivity n p h- 
In contrast to this, the low-T peak of n a i, is followed by a 
minimum around 80 K and a strong increase of K a b which 
eventually develops into a broad peak around 300 K that 
clearly exceeds the low-T peak. 

Several conventional mechanisms which may be in- 
voked to explain this unusual T-dependence and the re- 
sulting anisotropy can be dismissed. First, we can ex- 
clude effects of radiative heat transport since the opti- 
cal properties of La2Cu04 are almost isotropic in the 
relevant energy-range below hi/ < 0.1 eV [ij. Sec- 
ond, electronic contributions to the thermal current can 
be excluded since the material is insulating. Therefore, 
phononic heat transport is expected to dominate. How- 
ever, and third, a scenario solely based on thermal con- 
duction by acoustic phonons is very unlikely since an in- 
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FIG. 1: Thermal conductivity of La2Cu04 (bottom) and 
Lai.8Euo.2Cu04 (top). Full circles: k c . Open circles: K a b- 
Open squares: K mag . Solid line: fit according to Eq. Q 
Dashed line: K a b, P h- 

crease of scattering usually causes a negative slope of 
k(T) for intermediate and higher T, in contrast to our 
observation. Hence, either an unusual scattering process 
acts on the acoustic phonons or an additional contribu- 
tion apart from the conventional phononic background 
must be present for K a b- 

We now turn to the upper panel of Fig. ^ where K a b{T) 
and k c (T) of Lai.gEuo.2Cu04 are shown. Both curves 
strongly resemble our findings for La2Cu04, yet exhibit- 
ing obvious differences: the low T-peaks of K a b and k c 
are slightly larger and more sharply shaped than in the 
case of La2CuC>4. Furthermore, a step-like anomaly is 
present at Tlt ~ 135 K. Above Tlt k c remains almost 
constant and stays below the value for the undoped case, 
while K a b also exhibits a high-T maximum at T w 270 K 
which is even larger than in La2Cu04. 

The difference between k c of Eu-doped and of pure 
La2Cu04 can be attributed to a difference in phononic 
heat conduction. Upon doping La2Cu04 with Eu, en- 
hanced scattering of phonons reduces k c for T > Tlt, 
where both compounds have the same structure. The 
anomaly at Tlt (Lai.sEuo^CuCU) signals the transition 
to a new structural phase for T < Tlt where n p h is en- 
hanced. The structural peculiarities of rare earth doped 
cuprates are well known |l4j . Their influence on K p h will 
be discussed in detail in a forthcoming paper 0] . 

Despite the increase of scattering of phonons due to 
Eu-impurities in Lai.sEuo.2Cu04 the high-temperature 
peak of K a b is of larger magnitude than in La2Cu04. This 
excludes the peak to stem from the previously mentioned 



heat transport by acoustic phonons. The high-T peaks in 
K a b of both compounds therefore must originate from a 
different heat transport channel. In principle, two differ- 
ent kinds of excitation could generate such an additional 
heat current. One possibility is heat transport by dis- 
persive optical phonons, the other one is thermal current 
carried by magnetic excitations. Magnetic heat trans- 
port is the most likely explanation, since the magnetism 
is truly two-dimensional and the refore could explain the 
observed anisotropy. This is not true for dispersive opti- 
cal phonons, which certainly are present in this material, 
but are found along all crystal axes . 

This result is corroborated by the strong suppression of 
the high-T peak in hole-doped La? Y Sr Y CuQ4 0], even at 
a very low Sr-content of x 0.01 [l(j- On the one hand 
such an effect of doping is not to be expected in the case 
of a high-T peak originating from optical phonons, since 
the lattice impurities induced by the Sr-ions are similar 
to the Eu-impurities discussed above. On the other hand 
the strong frustration of antiferromagnetism upon dop- 
ing of mobile holes (cf. Ref. 0] and references therein) 
would provide for a satisfying explanation for the sup- 
pression of a peak of magnetic origin. Note in this context 
that Eu-doping leaves the Cu02-planes and therefore the 
magnetism almost unaffected. Hence we conclude, that 
the high-T peak of K a b originates from magnetic excita- 
tions which propagate only within the Cu02-planes. n a b 
therefore consists of a usual phonon background and a 
magnon contribution K mag while k c is purely phononic. 

Upon applying a magnetic field of 8 T no significant 
changes of k were detected. This is, however, consistent 
with a magnetic origin of the high-T peak since the corre- 
sponding Zccman-energy is orders of magnitude smaller 
than the magnetic exchange coupling J. 

In order to extract K mag from K a b we make use of the 
anisotropy of k assuming that the phononic part K a b,ph 
of K a b is roughly proportional to k c which is justified 
by only weakly anisotropic elastic constants |l5|. Since 
the magnetic contributions roughly follow a T 2 -law (see 
below) they are expected to be negligible in the range 
of the low-T peak. Therefore, for La2Cu04 as well as 
for Lai.gEuo.2Cu04 a reasonable estimate of ft a fc,ph is 
achieved by scaling the corresponding data for k c such 
as to match its low-T peaks with that of K a b- In Fig. 
the data thus obtained for K a f>,ph of both compounds are 
represented by dashed lines. K mag was extracted by sub- 
tracting K a b, P h from K a b (open squares in Fig. QJ. This 
procedure involves significant uncertainties only at low T 
where K mag <C K a b.ph- At intermediate and higher T the 
relative uncertainties are smaller and rather concern the 
magnitude of K a b,ph. 

than its slope. Therefore we assume 
that in this T-range the difference between the extracted 
and the true K mag is roughly constant. 

Following standard linearized Boltzmann-theory we 
have x l = \j^yz J ^kCk^f^k^k for the 2D thermal 
conductivity of a single magnon dispersion branch (la- 
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Sample Z mag (A) Fit interval (K) 



S La 2 Cu0 4 


558 ± 140 


70- 


158 


S Tifli a Will n oCllOi 


1157 ± 60 
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P La2Cuo.995Zno.00504 


515 


70 


158 


P La2Cuo.992Zno.00804 


435 


82 


164 


p La2cuo.99zno.0104 


266 


100 


-173 


p La2cuo.98zno.0204 


153 


118 


-200 


p La2cuo.95zno.0504 


49 


141 


-223 



TABLE I: Fit parameter of ft mag for single (S) and poly- 
crystals (P). The relative errors for Z mag of the polycrystals 
amount to about 10%. In all cases the errors have been ob- 
tained using different estimates for re p h. 




100 200 300 100 200 300 
Temperature (K) Temperature (K) 



belled by i), where v^, l^, £k and rik denote velocity, 
mean free path, energy and Bose- function of a magnon. 
Note that of a three-dimensional ensemble of planes, 
as realized in La2CuC>4, is given by K z mag = |>f*, where 
c = 13.2 A is the lattice constant of La2CuC>4 perpen- 
dicular to the planes. Then the total re mag results from 
summing up n l mag of each magnon branch. 

In order to calculate n l mag we approximate the magnon 
dispersion relation of the two branches i = 1 , 2 with 
the 2D-isotropic expression ej. = tk — yAf + i^ v ok) 2 , 
which describes the dispersion observed experimentally 
[rH llif for small values of k. Here, vq is the spin wave 
velocity while Ai and A2 denote the spin gaps of each 
magnon branch. For clarity, we define the characteristic 
temperature 9m = (?ivoy/ir) / (ciks) where a = 3.8 A is 
the lattice constant of the Cu02-planes 0] . Assuming a 
momentum independent mean free path, i.e. Ik = (mag, 
we obtain for each magnon branch: 

i _ vok B W T 2 [ 2 1 dn{x) 

mag " 2a 2 c G 2 M J V °' 5 dx { ) 

2:0,1 

Here, the integral is dimensionless but temperature de- 
pendent via xq,{ — Ai/(fcsT) and x maK . The upper 
boundary ic max , however, may be set to infinity with- 
out affecting the fit at temperatures T < 300 K. Since 
from neutron scattering experiments vo « 1.287- 10 5 m/s 
as well as A 1 /k B ~ 26 K and A 2 /fc B w 58 K [l3 are 
well known quantities, Z mag is the only unknown param- 
eter in Eq. ^ Assuming Z mag to be T- independent at low 
T, we can use Eq. order to check the T-dependence 
of our experimental ft mag , which for T > A1.2 should 
roughly be K mag oc T 2 , as the integral is only weakly 
T-dependent in this T-range. Furthermore, Z mag can be 
extracted from the data. For T > 250 K interactions be- 
tween the magnons become important j20| which requires 
a renormalization of the spin wave parameters and there- 
fore we restrict the application of Eq. ^ with constant 
l mas to T < 250 K. 



FIG. 2: Open circles: Thermal conductivity k of 
La 2 Cui_ z Zn z 4 polycrystals (z = 0, 0.005, 0.01, 0.05) as 
a function of T. Solid lines: extrapolated K p h- 



Fitting Eq. to the experimental data we allow for 
an additive shift of the K mag -curve which yields a fur- 
ther free parameter apart from Z mag and accounts for the 
aforementioned uncertainties in the magnitude of « mag 
(21) . For both compounds satisfactory fits were obtained 
at intermediate T-ranges. These are listed in Table 
together with the resulting values for Z mag . The fitting 
curves are represented by solid lines in Fig. ^ While 
the slight deviations between the fitted and experimen- 
tal data towards low T are due to the uncertainties in 
ftab, ph hi this range, the deviations at high T can be un- 
derstood in terms of the T-dependence of l mag due to 
enhanced magnon-magnon scattering. 

The data are consistent with a constant l mag for T 
both, within the fit interval and below, indicating that in 
this range T-dependent scattering processes like magnon- 
magnon scattering or even magnon-phonon scattering 
may be discarded. Therefore, relevant processes seem 
to be sample-boundary scattering or scattering at static 
magnetic defects. For Lai.sEuo.2Cu04 and La2Cu04 we 
find i mag « 1160 A and l ma g ~ 560 A, respectively. 
Since these values are far too small to correspond to the 
crystal dimensions, magnon-defect scattering is the most 
likely candidate. This conclusion is consistent with the 
fact that K mag and i mag are quantitatively different for 
Lai 8 Euo.2Cu04 and La 2 Cu04: since the magnetic prop- 
erties of both compounds are expected to be identical in 
essence, unequal K mag can arise only due to a difference 
in densities of the magnetic defects that restrict l mag . 

In order to check our analysis quantitatively it would 
be desirable to measure the density of static magnetic 
defects of our crystals. In lack of such methods we have 
performed measurements of k on samples with a well de- 
fined density of magnetic defects. Such defects can be 
induced in La2CuC>4 by substituting a small amount of 
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fore, this result confirms the above quantitative analysis 
of Kmag of the single crystals based on Eq. ^ 

In conclusion, both, qualitatively as well as quantita- 
tively our results strongly suggest, that the high-T peak 
observed in K a \, of doped and undoped La2Cu04 arises 
due to magnon heat transport which is confined to the 
CuC>2-planes. We note, that in general a large K mag is 
rarely observed at high T and requires a unique synergy 
of magnon-phonon coupling 26] , large Z mag and high spin 
wave velocity vq. Its realization in La2Cu04 may lead to 
future use of magnetic heat transport as a tool to study 
the interactions of magnetic excitations with other quasi- 
particles like holes or phonons. 

We acknowledge support by the DFG through SP1073. 



FIG. 3: Main panel: Magnon thermal conductivity «; mag of 
La 2 Cui_ z Zn z 4 (z = 0, 0.005, 0.008, 0.01, 0.02, 0.05) as a 
function of T (open circles). Solid lines: fits according to 
Eq. Inset: Z mag as a function of 1/z in units of lattice 
constants a. Solid line: fit line through origin. 



Cu 2+ -ions by non-magnetic Zn 2+ -ions. Representative 
results on our polycrystalline samples of La2Cui_ z Zn z 04 
are presented in Fig. [3 As is expected for doping of 
static structural and magnetic defects, the Zn-impurities 
lead to a gradual suppression of both, the phonon as well 
as the magnon contribution to k |25|. Due to the poly- 
crystalline nature of our samples anisotropic information 
on k is averaged over. We therefore estimate the phonon 



contributions Kp^ y by fitting k right of its maximum by 



Uph = a/T + (3 and by extrapolating this fit towards high 
T (solid lines in Fig.EJ. In turn, K%£^ on the polycrystals 
is obtained by Kjifag — K ~ K ph- 

Note, that the measured is smaller than the in- 

trinsic K mag of these compounds by the factor of 2/3 due 
to averaging over all three components of the k tensor 
[13. As shown in Fig. K mag = (3/2)«;P ol | systemat- 
ically decreases with increasing Zn-content as it is ex- 
pected for static defects. Analyzing these data we pro- 
ceed analogous to the undoped case by using Eq. ^ How- 
ever, we neglect a slight reduction of the spin wave ve- 
locity [2^ and changes of the spin gaps induced by the 
Zn-ions. These effects lead to corrections smaller than 
the experimental error. The fits are represented by the 
solid lines in Fig. |3 the resulting Z mag and the fit in- 
tervals are reproduced in Tabled For typical transport 
experiments it is expected that the mean free path of 
the quasi-particles is proportional to the reciprocal de- 
fect concentration, i.e., to 1/z. We therefore plot Z mag 
as a function of 1/z in the inset of Fig. [3] The linear 
scaling between Z mag and 1/z is clearly confirmed. From 
a linear fit including the origin (solid line) we find that 
( mag « 0.74a/ z. Since 1/z is equal to the mean unidirec- 
tional distance of Zn-ions, i.e. static magnetic defects, 
( mag gives a direct measure of these distances [24j . There- 
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